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INTRODUCTION 


The  first  Continuous  Automatic  Single  Base  Propellant  production 
Line  (Project  CASBL)  will  go  on  stream  at  Radford  Army  Ammunition  Plant. 
One  of  the  critical  problems  which  must  be  solved  to  provide  for  the 
operation  of  this  propellant  production  line  is  the  development  of  a 
rapid  and  accurate  method  of  charge  weight  estimation  to  allow  on-line 
loading  of  gun  charge  packages  without  the  delays  encountered  in  using 
proving  ground  gun  firings  to  establish  propellant  charge. 

An  on-line  gun  interior  ballistic  computer  model  has  been  proposed 
to  replace  the  gun  in  the  charge  establishment  role.  The  placement  of 
such  a model  in  a continuous  propellant  manufacturing  plant  is  illustra- 
ted in  Figure  1.  Figure  1 is  a flow  diagram  showing  the  relationships 
between  the  propellant  line,  the  propellant  line  process  controller, 
and  the  interior  ballistic  model.  Raw  material  goes  in  at  the  top  of 
the  line;  the  fabricated  propellant  grains  come  out  the  bottom,  are 
packed  into  gun  charges,  and  sent  to  storage.  Samples  of  the  production 
are  taken  off  at  various  parts  of  the  line,  subjected  to  chemical  and 
physical  tests,  and  the  information  from  these  tests  is  sent  to  a 
process  controller  which  in  turn  controls  the  operation  of  the  machines 
on  the  line  in  order  to  produce  a product  of  acceptable  quality.  Prior 
to  charge  packing,  propellant  samples  are  fired  in  the  closed  chamber. 
Information  from  the  closed  chamber  tests  and  from  the  chemical  and 
physical  tests  go  to  the  gun  interior  ballistic  model;  information  from 
the  gun  model  then  goes  to  the  process  controller  and  the  charge 
packing  station. 

Functional  details  of  the  gun  model  are  illustrated  in  Figure  2. 
Propellant  chemical  composition  data  from  chemical  analysis  will  be  pro- 
cessed by  a gun  thermodynamic  model  such  as  the  BLAKE  code*  to  produce 
propellant  thermodynamic  data  for  use  by  the  gun  simulation  model.  Pro- 
pellant combustion  characteristics  and  energy  data  will  be  provided  to 
the  gun  simulation  model  by  the  closed  chamber  data  reduction  program. 
Physical  examination  of  the  propellant  will  provide  the  gun  model  with 
propellant  density  and  grain  dimensions.  Gun  dimensional  data  and  pro- 
jectile weight  will  also  be  provided  to  the  gun  model.  The  gun  simulation 
model,  in  turn,  will  provide  muzzle  velocity,  maximum  breech  pressure, 
and  other  interior  ballistic  trajectory  data.  These  gun  performance  data 
will  then  be  compared  with  the  propellant  acceptance  criteria  and  a decis- 
ion made  to  accept  or  reject  the  propellant  lot.  If  the  propellant  lot 
is  rejected,  this  information  is  sent  to  the  process  controller  for 
appropriate  action.  If  the  propellant  lot  is  accepted,  a charge  weight 
will  be  computed  and  this  information  sent  to  the  process  controller. 


Eli  Freedman,  nA  Brief  Users  Guide  for  the  BLAKE  Program, " Ballistic 
Research  Laboratories  Interim  Memorandum  Report  249,  July  1974 . (Interim 
report  no  longer  available) 
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Fig.  l.Flow  Diagram  for  Continuous  Propellant  Manufacturing  Plant 
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Flow  Diagram,  Interior  Ballistic  Model  for  Continuous  Propellant  Testing 


The  ballistic  acceptance  plans  for  the  155-mm  howitzer  and  for  the 
175-mm  gun  charges  are  given  in  Tables  I and  II.  For  the  175-mm  gun, 
this  information  was  obtained  from  Reference  2 and  for  the  155-mm 
howitzer,  from  Reference  3.  It  will  be  noted  that  the  maximum  standard 
deviation  in  muzzle  velocity  mentioned  in  these  tables  is  the  round-to- 
round  standard  deviation  produced  by  measurement  errors,  variations  in 
projectile  and  gun  dimensions,  and  variations  in  the  loaded  charges. 

Table  I.  Ballistic  Acceptance  Plan  for 
175-mm  Gun  Propellant  Charge,  M86A2 

Sample  - 5 rounds  @ 70°F  (2l°C) 

5 rounds  @ 140°F  (60°C) 


70°F  (21°C) 

Velocity  Limits:  2940-3060  f/s  (896-933  m/s) 

Maximum  Standard  Deviation:  10  f/s  (3  m/s)  0.33% 

3a:  30  f/s  (9  m/s)  0.99% 

140°F  (60°C) 

Velocity  Limits:  3015-3130  f/s  (919-954  m/s) 

Maximum  Standard  Deviation:  16  f/s  (5  m/s) 

Maximum  Zone  3 Pressure:  48,400  psi  (334  MPa) 

Maximum  Excess  Charge  Pressure:  58,000  psi  (400  MPa) 


Table  II.  Ballistic  Acceptance  Plan  for 
155-mm  Howitzer  Propellant  Charge,  M4A2 

Sample  - 5 rounds  @ 70°F  (2l°C) 

5 rounds  @ 140°F  (60°C) 

Velocity  Limits:  1825-1875  f/s  (556-572  m/s) 

Maximum  Standard  Deviation:  9.2  f/s  (2.8  m/s)  0.50% 

3a:  27.6  f/s  (8.4  m/s)  1.50% 

Permissible  Individual  Pressure  @ 70°F  (21°C):  40,500  psi  (.279  MPa) 
Permissible  Individual  Pressure  @ 140°F  (60oC}:  44,100  psi  (304  MPa) 


2 

Fred  J.  Fitzsimmons,  "Concept  of  Scope  of  Work  for  PEMA  Project 
5774186  Acceptance  of  Propellant  Via  the  Continuous  Process  ( PROJECT 
AUTOCAP) Product  Assurance  Directorate,  Picatinny  Arsenal  Report  No. 
SMUPA-QA-A-P-55-73 , December  1972. 

j 

"Proving  Ground  Acceptance  Test  Procedure  for  Artillery  Propellants , 
Charge,  Propelling,  155-mm  M4A2,"  Aberdeen  Proving  Ground  Supplement 
P-155H- 7J, Cl,  13  January  1975  and  5 August  1974. 
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DESCRIPTION  OF  INTERIOR  BALLISTIC  METHODS 


Two  interior  ballistic  methods  were  used  in  this  study:  the  differ 
ential  coefficient  method  and  the  ballistic  similitude  method.  In  the 
differential  coefficient  method,  the  muzzle  velocity  and  the  maximum 
breech  pressure  are  predicted  by  two  equations: 


AV 

V 


al  + bl 


m + c (^fi)  + d 

RQ  1 1 RQj  ] 


ARF 

RF 


,ARF 

+ ei 


+ f. 


ACW 

cw 


ACW. ' 

1 cw- 


Cl) 


AP 

P 


a„  + b 


m + c cm,  . d ml 

2 RQ  2 ^ RQJ  2 RF 


.ARF.  r 

+ e2  RF  + f 


ACW 
2 CW 


.ACW, ' 
g2  1 CWJ 


(2) 


where:  AV 

V 

ARQ 

RQ 

ARF 


predicted  muzzle  velocity  minus  standard  muzzle  velocity 
standard  muzzle  velocity 

experimental  relative  quickness  minus  standard  relative 
quickness 

standard  relative  quickness 

experimental  relative  force  minus  standard  relative  force 


RF  = standard  relative  force 

CW  = standard  propellant  charge  weight 


ACW  = experimental  propellant  charge  weight  minus  standard 
propellant  charge  weight 

AP  = predicted  maximum  breech  pressure  minus  standard  maximum 
breech  pressure 


P = standard  maximum  breech  pressure 


a^,a^  = constant  terms 

b^,b^,etc  = differential  coefficients 

For  small  perturbations  in  relative  quickness,  relative  force  and 
charge  weight,  the  constant  and  second  order  terms  should  disappear.  In 
fact,  good  estimates  for  these  differential  sensitivity  coefficients  may 
be  made  from  a perturbation  analysis  of  either  the  ideal  gun  energy  con- 
servation equation  or  from  a suitably  more  complex  interior  ballistic 
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model.  In  practice,  more  accurate  results  are  usually  obtained  from  an 
empirical  correlation  of  the  experimental  data.  Higher  order  terms  in 
Equations  1 and  2 are  justified  for  those  cases  where  the  perturbations 
are  no  longer  small. 

The  differential  coefficient  method  has  been  used  previously  for 
charge  assessment.  Results  of  the  application  of  the  method  have  been 
reported  by  Jackson*^  and  by  Auerbach^. 

The  ballistic  similitude  method  was  originally  developed  by  Bennett^ 
in  1921  for  the  purpose  of  constructing  gun  interior  ballistic  tables. 
Since  in  those  precomputer  days  all  differential  equations  had  to  be 
integrated  by  hand,  simplification  of  the  governing  equations  was 
essential.  In  this  method,  two  unknown  parameters  (ballistic  quickness 
and  velocity  factor]  are  determined  by  simulation  of  a known  firing  and 
a convenient  change  of  the  independent  variable  in  the  derived  equations 
makes  it  possible  to  evaluate  these  unknowns  independently.  Hitchcock"7 
in  1956,  integrated  the  differential  equations  on  a digital  computer  to 
generate  an  extended  set  of  interior  ballistic  tables  to  meet  the  need 
to  deal  with  modern  weapons.  In  the  early  1960's,  Frankie  revised  this 
method  to  handle  weapons  with  two  propellants  (one  single-perforated 
and  one  seven-perforated  grain] . It  was  then  programmed  for  the  digital 
computer  so  that  instead  of  generating  tables,  each  interior  ballistic 
trajectory  problem  would  be  solved  as  an  individual  case.  Experimental 
maximum  breech  pressure  and  muzzle  velocity  are  needed  by  the  method  of 
compute  the  ballistic  quickness  and  velocity  factor  to  be  used  in  subse- 
quent new  trajectory  calculations.  The  ballistic  quickness  factor  is 
defined  by  the  following  equation: 


q 


q^CV)0'4 5 * 7  (p')°-5  (RF)£1  CR<ne2  Cl  - Kq(l  - | 

s 

WA  Cl  + e) 


4 

Wendell  F.  Jackson,  "Uniformity  of  Dupont  "NH"  Powder  for  the  155-mm 
Gun,"  Dupont  R eport  B L-204-7-6,  Burnside  laboratory,  November  1939. 

5 E.  Auerbach,  "The  Use  of  Closed  Bbmb  Data  for  Estimating  Charge  Weights 
for  105-mm  Howitzer,"  Picatinny  Arsenal  R eport  No.  1881,  May  1952. 

Q 

Albert  A.  Bennett,  "Tables  for  Interior  Ballistics,"  Ordnance  Depart- 
ment Document  No.  2039,  Washington , DC,  April  1921. 

7 

Henry  P.  Hitchcock,  "Tables  for  Interior  Ballistics , " Ballistic 
Research  laboratories  R eport  No.  993,  September  1956.  (AD  #134080) 
and  "Notes  on  Tables  for  Interior  Ballistics,"  Ballistic  Research 
laboratories  Technical  Note  No.  1298,  February  i960.  (AD  #234734) 
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Where:  q = 

ballistic  quickness  factor 

*1  = 

empirical  factor  determined  from  a known  firing 

CV  = 

chamber  volume 

P'  = 

reduced  projectile  weight 

e = 

weight  ratio 

A 

bore  area 

W 

propellant  web 

E1  = 

relative  force  exponent 

e2  = 

relative  quickness  exponent 

RQ"  = 

g 

web-independent  relative  quickness  after  Krier  and  Shimpi 

Kq  = 

q adjustment  factor  to  improve  charge  weight  correction 
(after  Heppner^} 

c = 

propellant  weight 

c = 

s 

propellant  weight  used  in  calibration  firing. 

The  ballistic  velocity  factor  is  defined  by  the  following  equation: 


r = 

rv  0.5  0.5 

rQ  (RF)  (1  - Kr(l  - f ))  (4) 

^ s 

Where:  r = 

ballistic  velocity  factor 

r = 

0 

empirical  factor  determined  from  a known  firing 

g 

gravitational  constant 

Kr  = 

r adjustment  factor  to  improve  charge  weight  correction 
(after  Heppner^) 

R.  Krier , M.  Means,  and  S.  Shimpi,  "Closed  Bomb  Data  Used  to  Predict 
Interior  Ballistics ,"  University  of  Illinois  Interim  R esearch 
Memorandum,  March  1974. 

Leo  D.  Heppner,  "Final  Report  of  Special  Study  of  an  Electronic 
Computer  Program  for  Interior  Ballistics , " Test  and  Evaluation  Cbmmand 
Report  No.  DPS-1.711,  /berdeen  Proving  Ground,  Maryland,  July  1965. 
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PROCEDURE  FOR  THE  ASSESSMENT  OF  THE  DIFFERENTIAL  COEFFICIENT  METHOD 


The  procedure  used  to  determine  the  differential  coefficients  and 
to  assess  the  method  was  as  follows: 

1.  Propellant  lot  acceptance  firing  data  from  firings  conducted 
at  Aberdeen  Proving  Ground  and  Jefferson  Proving  Ground  were  obtained 
and  selected  portions  of  these  data  were  stored  on  magnetic  tape 
cassettes  for  processing  by  a small  minicomputer. 

2.  Propellant  lot  description  data  were  obtained  from  Radford 
Army  Ammunition  Plant.  Selected  portions  of  these  data  were  also  stored 
on  magnetic  tape  cassettes. 

3.  The  coefficients  for  the  differential  coefficient  method. 
Equations  (1)  and  (2) , were  obtained  from  the  data  for  a selected  range 
of  propellant  lots  using  a least-squares  regression  technique.  For  the 
least-squares  analysis  and  statistical  tests,  standard  procedures  were 
followed10’  . 

4.  Once  the  coefficients  were  obtained  for  the  equations,  the 
muzzle  velocities,  maximum  breech  pressures,  and  recommended  charge 
weights  were  computed  either  for  the  range  of  lots  used  in  determining 
the  coefficients  or  other  ranges  of  lots  used  in  the  determination. 

5.  The  differences  between  the  predicted  and  measured  values  were 
computed  for  each  lot.  For  a range  of  propellant  lots,  the  values  of 
the  average  error  in  prediction,  the  standard  deviation  of  the  error, 
three  times  the  standard  deviation  of  the  error,  and  the  prediction 
confidence  intervals  were  computed.  The  prediction  confidence  interval 
brackets  the  predicted  value  such  that  one  is  99  percent  (as  computed 
here)  confident  that  the  next  observed  value  will  fall  within  the 
interval . 

Various  forms  of  the  differential  coefficient  equations  were  used 
in  this  study:  (1)  first  order  terms  with  no  constant  term,  (2)  first 

order  terms  with  constant  term,  (3)  first  and  second  order  terms  with  no 
constant  term,  and  (4)  first  and  second  order  terms  with  constant  term. 

The  155-mm  howitzer  and  175-mm  gun  propellant  lot  ranges  used  in 
this  study  are  listed  in  Table  III.  In  the  first  column  of  the  table, 
a symbol  designates  the  range  of  lots  used  in  the  study.  The  second 


", Experimental  Statistics  - Section  1,"  Ordnance  Engineering  Design 
Handbook,  ORDP-20-110,  June  1962,  Chapters  5 and  6 . 

H R.  L.  Anderson  and  T.  A.  Bancroft,  "Statistical  Theory  in  Research, " 
McGrow  Hill,  New  York,  1952 , Chapter  13. 
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column  designates  the  gun.  The  third  column  designates  the  range  of 
lots.  The  fourth  column  shows  the  reference  lot;  the  fifth  column,  the 
number  of  lots  in  a designated  lot  range;  and  the  sixth  column,  the 
number  of  maximum  breech  pressure,  muzzle  velocity  data  sets  used  in  the 
analysis.  This  last  number  is  greater  than  the  number  of  lots  for 
firings  because  more  than  one  charge  weight  per  lot  was  used. 

Lot  Range  C was  the  experimental  set  of  175-mm  propellant  lots 
whose  ballistic  properties  spanned  the  range  of  properties  expected 
during  propellant  production.  These  lots  were  used  to  predict  the 
performance  of  the  production  lots.  Results  from  Lot  Range  C predictions 
have  been  reported  previously  . Also,  six  lots  were  picked  from  Lot 
Range  A such  that  the  relative  force  and  relative  quickness  values 
exhibited  the  same  range  of  values  as  Lot  Range  A.  This  limited  set  of 
lots  is  designated  Lot  Range  AA.  The  same  was  done  for  Lot  Range  B to 
produce  Lot  Range  BB.  The  same  was  also  done  for  155-mm  Lot  Range  X to 
produce  Lot  Range  XX  except  that  only  five  lots  were  used.  Lot  Range  SX 
was  a set  of  experimental  155-mm  pn  pell ant  lots  whose  ballistic  properties 
spanned  the  range  of  properties  expected  during  propellant  production. 

These  lots  were  also  used  to  predict  the  performance  of  the  production 
lots . 


Selected  samples  of  the  better  results  are  listed  in  Table  IV.  In 
this  table  are  listed  the  correlation  lot  range,  that  is,  the  range  of 
lots  for  which  the  set  of  differential  coefficients  were  calculated; 
the  range  of  lots  for  which  muzzle  velocities,  maximum  breech  pressures, 
and  recommended  charge  weights  were  predicted;  the  gun;  and  the  parameter 
being  predicted.  This  is  followed  by  the  coefficients,  the  average  error, 
the  standard  deviation,  three  times  the  standard  deviation,  and  the 
99-percent  confidence  interval.  The  last  four  items  are  expressed 
as  percentages  of  the  standard  value.  Table  V gives  these  standard 
values  for  each  of  the  three  reference  lots  used  in  the  study.  Plots  of 
the  predicted  muzzle  velocity  or  maximum  breech  pressure  vs  the  measured 
muzzle  velocity  or  maximum  breech  pressure  for  some  of  the  differential 
coefficient  predictions  are  illustrated  in  Figure  3 through  Figure  6 for 
the  175-mm  gun  propellant  lots  and  Figure  7 through  Figure  10  for  the 
155-mm  howitzer  propellant  lots.  In  these  figures,  the  straight  line 
represents  the  locus  of  points  on  which  the  plotted  values  would  fall 
if  the  predicted  values  were  equal  to  the  measured  values.  The  two 
clusters  of  points  represent  data  for  upper  zone  charge  firings  and 
excess  charge  weight  firings.  Also  on  these  plots  are  illustrated  the 
mean  error,  the  standard  deviation  of  the  error,  and  the  differential 
coefficient  equation  used  to  predict  the  values  of  either  muzzle 
velocity  or  maximum  breech  pressure. 


12 

Paul  G.  Baevy  Ingo  W.  May , and  Jerome  M.  Frankie > "A  Comparison  of 
Several  Predictive  Approaches  in  Charge  Establishment  for  Large 
Caliber  Artillery  Systems 3ft  Proceedings  of  11th  JANNAF  Combustion 
^f^^bumeni  §b  ^asa^ena^  California , September  1974 . 
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Table 

III.  Propellant  Data  Used 

in  Regression  Study 

Reference 

Number 

Number  of 

Symbol 

Gun 

Propellant  Lots 

Lot 

of  Lots 

P,V,  Data  Set: 

A 

175-mm 

RAD  65009  to  RAD  65574 

BAJ 

65512 

87 

181 

B 

175-mm 

RAD  65751  to  RAD  67994 

RAD 

65306 

122 

235 

c 

175-mm 

RAD  A-PEI-441-1  to 
RAD  H-PEI-441-1 

RAD 

65306* 

67994** 

8 

16 

AA 

175-mm 

RAD  65017,  65026, 
65029 

65215,  65414,  65528 

BAJ 

65312 

6 

12 

BB 

175-mm 

RAD  65733,  65736 
65982,  65275,  67481 
67566 

RAD 

65306 

6 

12 

X 

155-mm 

RAD  69337  to  RAD  69444 

RAD 

68308 

52 

96' 

XX 

155-mm 

RAD  69398,  69421, 

RAD 

C 

15 

69422,  69439,  69443 

68308 

J 

SX 

155-mm 

RAD-PE-441-J  to 

RAD 

15 

30 

RAD-PE-441-W 

68308 

*Closed  chamber  reference  lot 

**Gun  firing  reference  lot 
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Table  IV.  Errors  in  Prediction  of  Muzzle  Velocities,  Maximum  Breech  Pressures,  and  Recommended  Charge  Weights 
for  the  155-mm  Howitzer  and  175-mm  Gun  Propellant  Lots  Using  the  Differential  Coefficient  Method 
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Standard  Value 


Table  V.  Standard  Values  for  155-mm  and  175-mm 
Reference  Propellant  Lots 


Reference 

Maximum 

Pressure 

Muzzle  Velocity 

Charge  Weight 

Lot  No. 

Gun 

psi 

MPa 

f/s 

m/s 

lb 

kg 

BAJ  65312 

175-mm 

51,100 

352.3 

3000 

914.4 

55.56 

25.20 

RAD  65306 

175-mm 

47,200 

325.4 

3000 

914.4 

57.24 

25.96 

RAD  68308 

155-mm 

36,500 

251.7 

1850 

563.9 

13.28 

6.02 

PROCEDURES  FOR  THE  ASSESSMENT  OF  THE  BALLISTIC 
SIMILITUDE  METHOD  (HITCHCOCK  CODE) 

The  procedures  used  to  determine  adjustable  parameters  for  the  bal- 
listic similitude  method  and  to  assess  this  method  were  as  follows: 

1.  The  propellant  lot  acceptance  firing  data  and  the  propellant 
lot  description  data  were  stored  in  data  files  on  the  UNIVAC  1108  com- 
puter^ magnetic  drum  and  disc  storage. 

2.  The  Hitchcock  code  was  run  on  the  UNIVAC  1108  computer,  varying 
the  two  empirical  factors  (q^  in  Equation  (3)  and  r^  in  Equation  (4))  in 
the  code  until  the  predicted  muzzle  velocities  and  maximum  breech  pressures 
agreed  with  the  muzzle  velocities  and  maximum  breech  pressures  obtained 
from  the  firings  of  the  reference  propellant  lots. 

3.  With  the  empirical  factors  held  constant,  the  Hitchcock  code 
was  used  to  predict  the  muzzle  velocities  and  maximum  breech  pressures 
for  a range  of  propellant  lots. 

4.  The  differences  between  the  predicted  and  measured  values  were 
obtained  for  the  test  propellant  lots  and  then  the  average  error,  the 
standard  deviation  of  the  error,  and  three  times  the  standard  deviation 
of  the  error  were  obtained.  Confidence  limits  could  not  be  computed 
using  the  values  from  the  Hitchcock  code  because  the  method  of  obtaining 
confidence  limits  does  not  apply  to  non-linear  regression  models  such  as 
the  Hitchcock  code. 

In  Table  VI,  the  errors  in  predicting  muzzle  velocities  and  maximum 
breech  pressures  for  the  175-mm  gun  and  155-mm  howitzer  production  lots 
using  the  Hitchcock  code  are  listed. 
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MEASURED  VELOCITY  - f/s 

Fig.  5.  Prediction  of  Muzzle  Velocities  for  175-mm  Gun  Propellant  Lots  Using  Differential 
Coefficients.  Lot  Range  A Predicted  from  Second  Order  Correlation  to  Lot  Range  B 
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42000  48000  50000  54000  58000 

MEASURED  PRESSURE psi 

Fig.  6.  Prediction  of  Maximum  Breech  Pressures  for  175-mm  Gun  Propellant  Lots  Using  Differential 
Coefficients.  Lot  Range  A Predicted  from  First  Order  Correlation  to  Lot  Range  B 
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Fig.  8.  Prediction  of  Maximum  Breech  Pressures  for  155-mm  Gun  Propellant  Lots  Using  Differential 
Coefficients.  Lot  Range  X Predicted  from  First  Order  Correlation  to  Lot  Range  X 
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0.  Prediction  of  Maximum  Breech  Pressures  for  155-mm  Gun  Propellant  Lots  Using  Differential 
Coefficients.  Lot  Range  X Predicted  from  First  Order  Correlation  to  Lot  Range  XX 


In  Table  VI  are  listed  propellant  lot  ranges  in  which  the  predic- 
tions were  made,  the  gun,  the  parameter  predicted  (either  muzzle  velo- 
city or  maximum  breech  pressure) , the  four  adjustable  parameters  in  the 
Hitchcock  code  which  were  varied  in  this  study,  the  average  error,  the 
standard  deviation  of  the  error,  and  three  times  the  standard  deviation 
of  the  error.  The  last  three  entries  are  expressed  as  percentages  of 
the  standard  values.  In  this  study,  three  production  propellant  lot 
ranges  were  predicted,  Lot  Ranges  A and  B for  the  175-mm  gun  and  Lot 
Range  X for  the  155-mm  howitzer.  The  normal  values  for  the  four  ad- 
justable param^l;ers  as  determined  by  previous  simulations  with  the 
Hitchcock  code  ^ are  given  in  the  first  line  of  Table  VI.  For  the  next 
six  sets  of  runs,  in  each  of  which  muzzle  velocity  and  maximum  breech 
pressure  are  predicted  simultaneously,  the  relative  force  exponent  was 
used  at  its  nominal  value  of  0.1667  or  zero.  A zero  value  for  the 
relative  force  exponent  would  imply  that  variations  in  relative  force 
would  have  no  effect  on  predicted  results.  Indeed  slightly  improved 
predicted  results  are  obtained  for  some  of  the  correlations  shown  in 
Table  VI.  Such  spread  in  relative  force  approaches  the  accuracy  of  the 
closed  bomb  measurements. 

Plots  of  the  predicted  muzzle  velocities  or  maximum  breech  pres- 
sures vs.  measured  muzzle  velocities  or  maximum  breech  pressures  for 
175-mm  gun  Lot  Range  B,  using  nominal  values  of  the  adjustable  para- 
meters except  that  the  relative  force  exponent  was  zero,  are  illustrated 
in  Figures  11  and  12. 


SUMMARY  OF  RESULTS 


All  of  the  results  using  either  the  differential  coefficient  method 
or  the  ballistic  similitude  method  to  date  for  the  155-mm  howitzer  and 
175-mm  gun  production  and  experimental  lots  are  summarized  in  Table  VII. 
The  results  are  given  in  terms  of  standard  deviation  of  differences  be- 
tween predicted  and  measured  results  and  in  terms  of  the  99-percent 
Confidence  Interval,  with  both  values  expressed  in  terms  of  percentages 
of  standard  values.  The  numbers  given  for  the  standard  deviation  and 
99-percent  Confidence  Interval  are  expressed  as  a range  of  values 
between  a minimum  and  a maximum.  Not  all  of  the  differential  coefficient 
correlations  summarized  in  Table  VII  were  detailed  in  Table  IV.  It  is 
possible  to  compute  the  recommended  charge  weights  using  the  ballistic 
similitude  method,  but  this  has  not  yet  been  done.  In  general,  the 
standard  deviation  of  error  values  for  muzzle  velocity  using  the 
differential  coefficient  method  is  about  one-half  that  of  the  ballistic 
similitude  method.  For  maximum  breech  pressures,  the  ballistic  similitude 
method  gives  standard  deviation  error  values  approximately  three  times 
those  given  by  the  differential  coefficient  method. 
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Prediction  of  Muzzle  Velocities  for  175-mm  Gun  Propellant  Lots  (Lot  Range  B) 

Using  the  Hitchcock  Code 
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Fig.  12.  Prediction  of  Maximum  Breech  Pressures  for  175-mm  Gun  Propellant  Lots  (Lot  Range  B) 

Using  the  Hitchcock  Code 


Table  VII.  Summary  of  Results 

Standard  Deviation  99%  Confidence  Interval 


Gun 

Model 

Variable 

% Standard  Value 

% Standard  Value 

175-mm 

DFC 

Velocity 

0.23-0.58 

0.69-  4.06 

175-mm 

DFC 

Pressure 

0.70-4.21 

4.64-29.46 

175-mm 

DFC 

Rec  Chg  Wt 

0.44-0.65 

1.23-  6.56 

155-mm 

DFC 

Velocity 

0.18-1.11 

0.67-  3.23 

155 -mm 

DFC 

Pressure 

0.94-9.40 

3.58-27.50 

155-mm 

DFC 

Rec  Chg  Wt 

0.41-0.46 

1.09-  1.63 

175-mm 

H 

Velocity 

0.49-0.73 

- 

175-mm 

H 

Pressure 

3.15-4.70 

- 

155-mm 

H 

Velocity 

0.51-0.61 

- 

155-mm 

H 

Pressure 

2.95-3.37 

_ 

DFC  - Differential  Coefficient  Method 
H - Ballistic  Similitude  Method  (Hitchcock  Code) 

Rec  Chg  Wt  - Recommended  Charge  Weight 

EFFECTS  OF  CLOSED  CHAMBER  MEASUREMENT  ERRORS 


An  attempt  was  made  to  assess  the  effects  that  errors  in  closed 
chamber  measurement  would  have  on  predicted  muzzle  velocity.  Unpub- 
lished data15  from  Radford  Army  Ammunition  Plant  gave  closed  chamber 
measurement  errors  for  the  175-mm  gun  experimental  propellant  lots  RAD- 
A-PEI-441-1  to  RAD-H-PEI-441-1 . These  data  are  summarized  in  Table  VIII 
for  the  particular  data  reduction  techniques  used  in  determining  the 
relative  quickness  and  relative  force  published  in  the  propellant  des- 
cription sheets.  The  average  values  were  used  in  the  regression  analy- 
sis for  the  differential  coefficient  method.  The  differential  coeffi- 
cient method  used  to  predict  muzzle  velocities  was  a first  order  fit 
with  constant  term  to  Lot  Range  C (RAD  A to  RAD  H) . This  equation  is: 

~ = 0.0009  + 0.299  ^ + 0.539  (5) 

Listed  in  Table  VIII  are  the  minimum,  maximum,  average,  and  stand- 
ard deviation  values  for  relative  quickness  and  relative  force.  Also 
given  are  the  predicted  sample-to-sample  standard  deviations  in  muzzle 
velocities  using  the  differential  coefficient  method.  Six  samples  from 
each  lot  were  fired  in  the  closed  chamber  to  give  the  values  listed. 

nz 

Letter , "Project  AUTOCAP  (RAAP  Project  PE-441),"  Radford  Army  Ammuni- 
tion Plant,  26  September  1974. 
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Table  VIII.  Effects  of  Closed  Chamber  Measurement  Errors 
on  Predictions  Using  First  Order  Fit 

175-mm  Experimental  Propellant  Lots  RAD  A to  RAD  H 
6 Closed  Chamber  Firings 
IBM  Filtered  Reduction  Technique 


RQ  RF  Vel  SD 


Lot 

Min 

Max 

Ave 

SD 

Min 

Max 

Ave 

SD 

f/s 

m/s 

A 

99.8 

104.0 

102.8 

1.57 

100.2 

102.1 

101.0 

0.70 

25.4 

7.7 

B 

109.6 

112.2 

110.8 

0.98 

101.2 

103.4 

102.0 

0.73 

20.6 

6.3 

C 

94.7 

99.0 

97.2 

1.54 

100.3 

101.1 

100.7 

0.56 

19.6 

6.0 

D 

103.2 

109.3 

106.5 

2.46 

101.7 

102.5 

102.0 

0.25 

26.1 

7.8 

E 

92.2 

94.8 

92.9 

0.93 

99.5 

100.4 

99.8 

0.36 

14.1 

4.3 

F 

100.3 

104.5 

102.8 

1.66 

100.9 

101.7 

101.2 

0.31 

20.0 

6.1 

G 

82.5 

85.0 

83.4 

1.01 

94.6 

95.9 

95.3 

0.52 

17.5 

5.3 

H 

83.2 

86.7 

84.9 

1.19 

94.6 

95.5 

95.1 

0.30 

15.5 

4.7 

SD  - Round-to-Round  Standard  Deviation 

Note  that  the  standard  deviations  in  predicted  muzzle  velocity  range 
from  14.1  f/s  (4.3  m/s)  to  26.1  f/s  (7.8  m/s).  This  represents  an  0.5- 
to  0.9-percent  variation  of  the  standard  muzzle  velocity.  Since,  accord- 
ing to  Table  I,  the  maximum  round-to-round  standard  deviation  is  0.33 
percent  of  the  standard  muzzle  velocity,  it  is  apparent  that  these 
muzzle  velocity  standard  deviations,  if  they  had  appeared  in  gun  firings, 
would  cause  the  lots  to  be  rejected.  Since  such  round-to-round  standard 
deviations  were  not  observed  in  the  gun  firings  for  these  lots  (round- 
to-round  standard  deviation  in  muzzle  velocity  ranged  from  0.13  percent 
to  0.33  percent  according  to  Reference  14),  it  is  apparent  that  these 
closed  chamber  standard  deviations  are  mainly  measurement  errors  and  do 
not  reflect  propellant  property  variations  within  each  propellant  lot. 
While  these  measurement  errors  will  not  cause  round-to-round  standard 
deviation  errors  in  gun  firings,  they  will  seriously  affect  the  ability 
of  a model  to  properly  predict  the  results  of  a gun  firing. 

CONCLUSIONS 


From  the  results  of  this  study,  it  can  be  concluded  that  the  wholly 
empirical  differential  coefficient  method  will  predict  values  of  muzzle 
velocities  and  maximum  breech  pressures  better  than  the  ballistic  simi- 
litude method.  The  drawback  of  the  differential  coefficient  method  is 


Id 

~ "Product  Improvement  Test  of  Charge  Propelling 175-mm > M86A2  CEval 
uation  of  Propellant  Variables  on  Ballistic  Performance) Firing 
Record  No.  P-82419 > Aberdeen  Proving  Ground 3 Maryland y 1 April  1974. 
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that  it  will  require  experimental  firing  of  a set  of  propellant  lots 
whose  ballistic  properties  will  span  the  complete  range  of  properties 
expected  during  production. 

The  ballistic  similitude  method  can  be  used  for  start-up  with  a new 
propellant  and/or  a new  gun  since  the  method  requires  a minimum  firing 
of  one  lot  to  calibrate  the  model.  More  research  and  better  input  data 
are  needed  in  order  to  completely  replace  the  differential  coefficient 
method  with  the  ballistic  similitude  method  or  other  interior  ballistic 
models  of  greater  complexity. 

Closed  chamber  measurement  errors  seriously  affect  the  ability  of 
the  models  to  predict  the  ballistic  performance  of  guns.  Improvements 
in  both  precision  and  accuracy  of  closed  bomb  burning  rate  and  impetus 
parameters  are  essential.  Coupled  with  these  improvements,  burning  rate 
parameters  different  from  the  conventional  relative  quickness  (e.g.  vivacity 
or  linear  burning  rate)  may  yield  more  accurate  predictions. 
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